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SECTION I. 


INTRODUCTION 


n 


This report summarizes the work during a 12-month study 
to establish and apply analytical techniques for development 
of a vector wind profile gust model for the Space Transporta- 
tion System OFT Operations and Trade Studies. 

Jones (Ref. 1) summarizes the fundamental issue involved 
in application of gust models for evaluation of vehicle struc- 
tural response. Basically, vehicle structures can be considered 
to be either highly damped or slightly damped. Highly damped 
modes are associated with vehicle control and handling qualities; 
large response of these modes can be generated by imposing a 
single discrete gust with an appropriate length scale. On the 
other hand, large responses for slightly damped modes could be 
the cumulative effect of a sequence or cluster of gusts with an 
appropriate spacing or phasing; such a sequence of gusts is best 
treated with a power spectrum im^del. 

The application of either model depends on which type of 
response mode is dominant. The inaccuracy of using one model 
and neglecting the other should be evaluated and taken into 
account. For example, when power spectra are used, analysis of 
the effect of the non-Gautsian structure of turbulence is 
necessary; when a discrete gust model is used, simple cluster- 
ing of discrete gusts should be studied for possible resonance 
effects. 

Jones concludes that analysis of large disturbances re- 
corded during passage through patches of turbulence are usually 
associated with a single gust that stands out above the general 
level of roughness. It is implied that these large loads can 
be simulated with a discrete gust model. 

From the above discussion, it is clear that a discrete 
gust model can be useful in studying the response of damped 
vehicle structures to turbulence. The results of this study 
will be suitable for modeling of discrete gusts in vertical 
wind profiles. 

The body of this report is composed of three sections 
(II through IV). Section II describes basic and derived pro- 
perties of the data analyzed in this study. Section III 
describes the calculation of wind gust data from Jimsphere 
wind profiles and the ested>lishment of probability distribu- 
tions of gusts in four wavelength bands. Conclusions are 
presented in Section IV. 


1 


1 



SECTION II. DATA 


This section describes basic and derived properties of 
the data analyzed in this study. Basic properties include 
the number, type, and location of the wind profiles; derived 
properties include an evaluation of the response of the 
measurement system at small wavelengths. 


A. Data Sample 


The data consist of 1800 Jimsphere profiles (150 per 
month) from Cape Kennedy, Florida (Ref. 2 ) . The data were 
obtained from a Space Shuttle Level II directive that speci- 
fies the demonstration of vehicle design validity using 
150 Jimsphere wind profiles representative of each month. 

Three months (February, April, and July) were chosen for analy- 
sis in this study. The February and July data are representa- 
tive of the seasonal extremes at Cape Kennedy; the April data 
are representative of the transition between the extremes. 

The nximber of soundings for each month for each year of the 
sampling period is illustrated in Figure 1. 


B. Wind Measuring System Amplitude Response 


wind profile data used in this study were obtained with 
the Jimsphere system. Since the small wavelength perturbations 
observed in these profiles are the subject of a detailed 
analysis, it is appropriate to specify the accuracy of the 
system for small wavelengths. A measure of the accuracy is 
the amplitude respo; se, G(i), which is equivalent to the 
ratio A(X)/A*(X); where A*(X) is the true amplitude of a per- 
turbation in the wind profile at wavelength, X, and A(X) is 
the amplitude measured with the Jimsphere system. The ampli- 
tude response of the Jimsphere system is limited by the size 
of the balloon (2-meter diameter) , the balloon ascent rate 
(4-5 m/sec), the accuracy of the balloon tracking system 
(FPS-16) , and the data smoothing technique. The balloon po- 
sitions, determined every 0.1 second, are smoothed to provide 
mean positions at each 25-meter interval of ascent* 

Differences in position between alternate 25-meter levels 
indicate the mean wind for the corresponding 50-meter layer, 
and are reported as the wind at the 25-meter level in the mid- 
dle of the 50-roeter layer. Thus, the basic data analyzed here 
are wind speeds and directions for 50-meter layers, overlapping 
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YEAR 


February 


Figure 1. Distribution of February » April, and July 
Jimsphere Soundings from 150 per Month 
Sample (Ref. 2 . 
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by 25 meters. Only when at least 25 meters intervene between 
two layers (i.e., winds reported for levels at least 75 meters 
apart) can two winds be considered independent observations 
(Ref. 3 ) . 


Expressionr for the amplitude response, G{X), of the Jim- 
sphere system to wind perturbation wavelengths that are small 
relative to the length of the wind profile have be'^n derived 
by Luers and Engler fRef. 4), 


G(X) = 




and by DeMandel and Krivo (Ref. 5), 

r4TTW\ /50 tt 


G(X) = 


. /4ttw\ /50tt\ 
sin(-y-) sin(-^-) 

200w(J.) 


( 1 ) 


( 2 ) 


where 


S = smoothing interval = 75m 
X - wavelength (m) 

w = Jimsphere balloon ascent rate (m/s) 

As illustrated in Fig. 2 , the Jimsphere system does not 
measure wavelengths less than 50 meters; for X=90m, the 
measured amplitude is one-half the true amplitude. 
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SECTION III. WIND ANALYSIS 


This section describes the calculation of wind gust data 
from Jimsphere wind profiles and the establishment of pro- 
bability distributions of gusts. Sample estimates of the 
required parameters of the theoretical distribution functions 
are calculated and the variability of the parameters as a 
function of altitude, season, and gust wavelength range is 
established. Empirical functions for estimation cJ gust sta- 
tistics for other wavelength ranges are derived. 


A. Digital Filters 


Vector wind gust statistics and models are based on data 
that have been obtained from filtered wind profiles. The fil- 
tering process provides profile data that contain perturba- 
tions within a range of wavelengths that is suitable for simu- 
lation studies of space vehicle ascent through the atmosphere. 
The design and application of these filters are described 
below. 


1. Filter Design . The design of the digital filters 
is based on the Martin-Graham cosine rolloff model described 
by Demandel and Krivo (Ref. 6 ) . A set of numerical smoothing 
weights is calculated for a low-pass filter from the equation 


h(nT) 


sin(2TTf.nT) + sin(2Trf nT) 
t c 

2TrnT f 1 - 4n^T^ (f ^ - f^) ^ 


(3) 


where the filter design parameters are 

T = altitude interval of wind profile data 

n = weight index (-N, -N+1, ..., -1, 0, 1, ..., 
N-1, N) 

N = (NW-D/2 

NW = number of weights 

f = cutoff frequency = the highest frequency with 
associated amplitude passed with unity gain 

f. = termination frequency = the lowest frequency 
with associated amplitude passed with ..ero 



The center weight (n = 0) is given by: 


h 


0 



(4) 


When the weights, h^, have been determined, they are normalized 
by applying the constraint 


N 

E = 1- 

N=-l 


(5) 


Only (N + 1) weights are calculated since h^ = Since 

the filter function is symmetrical, no phase shift is 
produced. 

The use of digital smoothing weights results in the loss 
of the first and last N data points of the original profile. 
Thus the filtered wind profile has an altitude range that is 
reduced by 2NT compared to the original profile. 

The effective response of the low-pass filter, given the 
design parameters listed under equation (3) is 

N 

G (f) = h. + 2 X) h cos(2TifnT). (6) 

^ ^ n=l ” 


As the number of weights (NW) is increased, the response 
of the filter improves. However, computation time increases 
as does the number of points lost (the first and last N data 
points) . In this study, NW was chosen to minimize data loss 
while maintaining a reasonably accurate filter response. 

2. Filter Application . Jimsphere wind profiles from 
the surface to 20 km in component form (zonal and meridional) 
were decomposed into eight data bases by the filtering pro- 
cess diagrammed in Figure 3. Four of the data bases consist 
of low-pass profiles that can be used in analyses of steady 
state and wind bias profiles. The other data bases consist 
of high-pass profiles defined here as residual profiles; 
these profiles consist of perturbations with relatively small 
wavelengths that are of interest in evaluations of vehicle 
response. Gusts that are derived from residual profiles are 
the subject of the detailed statistical analysis described 
in subsequent sections of this report. 
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Low'Pass Filter III 


_J 


Low'Pass Profile III 
(LPIII) 

X > 997 m 


J - LPIII “ 

Residual Profile III 
90 < X < 997 m 


Low'Pass Filter IV 


Low-Pass Profile IV 
(LPIV) 

X > 6000 ID 


J - LPIV « 

Residual Profile IV 
90 < X < 6000m 


^Nominal low wavelength limit of Jimsphere system 
’Defined in text 


Figure 3 . Filtering Process 
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The design parameters and weighting functions of four 
low-pass filters and the altitude range of the filtered pro- 
files used in this study are listed in Table 1 . 

The method of calculating high-pass profiles by subtrac- 
tion of the low-pass filtered profiles from the original 
Jimsphere profile is equivalent to the execution of a high- 
pass filter. The effective amplitude response of the four 
high-pass filters that are appropriate for description of the 
upper end of the wavelength range of the residual profiles is 
illustrated in Figure 4 . The nominal high wavelength limit 
for each set of residual profiles is the wavelength at which 
the amplitude response of the corresponding filter is .50. 

A set of u component residual profiles calculated from 
the Jimsphere profile of 28 July 1965 (2259 Z) at Cape Kennedy 
is illustrated in Figure 5 . 


B. Definition of Gust 


The definition of gust used in this study satisfies the 
objective to provide data that are suitable for a detailed 
statistical analysis of singularities and quasi-sinusoidal 
perturbations that are often observed in Jimsphere wind pro- 
files. A statistical model of these gusts so defined will be 
developed that will be useful for certain types of flight 
simulations of space vehicle ascent through the perturbed 
atmosphere. 

According to the conventional approach, a gust profile 
is calculated by applying a high-pass digital filter to a 
Jimsphere profile; all the magnitudes in the filtered profile 
are defined as gusts. In this study, these magnitudes are de- 
fined as residuals; the maximum positive or negative ' residual 
in the vicinity of a specified reference altitude is defined 
as a gust. A formal definition of gust is given below. 

Let u* represent the zonal wind component at a specified 
reference altitude, in a residual profile. The zonal gust 

is defined as the maximum value of u* in the vicinity of alti- 
tude with like sign to u' at H . The altitude interval 

associated with the gust is defined as the gust length, L, 
which is calculated by taking the altitude difference of the 
zero crossings on either side of the gust; i.e., 

L = - Hi (7) 
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Table 




Filter Design Parameters, Filter Weighting Functions 
of Four Filters Used for Calculation of Residual 
Profiles, and Altitude Range of Residual Profiles 


f liter D esig n Parame t e 


filter 

I 

II 

in 

IV 

F ilter Weights 


ho 4 2 t 


T(m) 

25 

250 

25 

250 






20 

5 

50 

10 


1 

0.116360050 

0.112681533 

0.102183235 

0.086369542 

0.067415386 

0.047750214 

0.029618173 

0.014711243 

0.003949048 

-0.002560992 

-0.005394941 

-0.005565475 

-0.004229394 

-0.002423366 

-0.000884042 

0.000021198 

0.000259004 

0.000022211 

-0.000405784 

-0.000771288 

-0.00092SS30 


1.000000001 


.00034 


.00435 

.00004 


.00080 

.00050 


.00150 

0 


.000342 

11 

III 

IV 

0.203331671 

0.050406609 

0.084765087 

0.182602840 

0.050170253 

0.083178582 

0.130080937 

0.049465762 

0.078561135 

0.068650095 

0.048306755 

0.071321355 

0.020649325 

0.046715542 

0.062084219 

-0.003649032 

0.044722562 

0.051615690 


0.042365613 

0.040733073 


0.039688904 

0.030213801 

1.000000000 

0.036741958 

0.020715102 


0.033578388 

0.012714788 


0.030254595 

0.026828417 

0.023357771 

0.006479712 


0. Cl 9899321 
0.016507215 
0.013231923 
0.010119200 
0.007209211 
0.004535825 
0.002126107 
0.000000000 

1.000000000 


-0.001829786 

-0.003357689 

-0.004585029 

-0.005519564 

-0.006174897 

-0.006569761 

-0.006727193 

-0.006673649 

-0.006438065 

-0.006050919 

-0.005543301 

-0.004946033 

-0.004288852 

-0.003599693 

-0.002904062 

-0.002224550 

-0.001580449 

-0.000987520 

-0.000457877 

.333639299-12 

0.000381141 

0.000683858 

0.000909429 

0.001061717 

0.001146732 

0.001172171 

0.001146932 

0.001080646 

0.000983218 

0.000864417 


0.999999999 


Altitude Hinge of 

Hesiduil 

Hrofilet 

Filter 

^■in 



(kn) 

(Ml) 

I 

O.S 

19.50 

11 

1.75 

18.25 

111 

1.25 

18.75 

IV 

3.75 

16.25 
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Figure 4 . Amplitude Response of Four Digital Filters 
Used for Calculation of Residual Profiles 


ALT. KM 


UAvri cnctm raw 

90-H?0f1 90-997M 

а. 

10 . 

18. 

17. 

16. 

15. 

!«#. 

13. 
la. 
n. 

10. 

0 . 

8 . 

7. 

б . . 

5. 


3. 

• • . ! 

0 [iiin] 0 

-e 0 e -^*-2 9 ^ H 



u. M/f U. M/5 


CAPt Kr^^ADY RESIDUAL PROriLCS 
U COri\U NT 

9D-?M70ri 


90-6000M 




Figure 5. Cape Kennedy Residual Profiles 
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The altitudes of the zero crossings, H 2 and are calculated 
by linear interpolation according to 


= H. 


25 


u 


*2 “j-1 u! - “j-1 

“1 = “k+1 


( 8 ) 


where 


«2 


altitude of the first zero crossing for 
the upward scan 


= 

last value of u' with the like sign of u* 
at Hq when scanning upward* 

u! 

D 

=r 

first value of u' with sign opposite to 
sign of u' at Hq when scanning downward 

“i-1 


altitude of 

«1 


altitude of the first zero crossing for 
the downward scan 

'^k+l 

= 

last value of u' with like sign to sign of 
u' at Hq when scanning downward 

“k 


first value of u' with sign opposite of u* 
at Hq when scanning downward 

«k+l 

= 

altitude of 

Similarly, 

substitution 

the meridional gust component, v', is defined 
of v' for u' above. In most instances, the 


zonal and meridional component gusts defined in this manner do 
not occur at the same altitude. This altitude difference, AH, 
is a measure of the phase difference between the components. 

A schematic definition of gust is given in Figure 6 . 


1 


I 

I 


^The indices j and k increase upward. 
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C. Gust Statistical Parameters 


The calculation of theoretical gamma distributions of 
component gust and associated gust length requires that esti- 
mates be made of the parameters of the distribution. Methods 
for estimation of the parameters are described in this sec- 
tion. Examples are presented which illustrate the impact of 
parameter estimation on the calculated theoretical distribu- 
tions. Comparisons with observed distributions are presented. 
The variability of the parameters as a function of altitude 
and filter cut-off frequency is also described. 

1. Parameter Estimation . The parameters of the gamma 
distribution are estimated from sample statistics. The scal- 
ing parameter, 6/ is calculated from an estimate of the shape 
parameter, y, according to 


6 = Y/x (10) 


The parameter, y, can be estimated according to the 
moments method (M) 


Y = (x/o)* 


( 11 ) 


where x and o are the mean and standard deviation of the data 
sample . 

Alternatively, an estimate of y can be obtained from 
either of two maximum likelihood procedures; according to 
Thom (Ref. 7), y is given by (MLI) 


Y 


^ (1 + ✓ 1 + 4A/3 ) 


(12) 


where A ■ ln(x) - In^x) 


(13) 
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According to Bury (Ref. 8), Y is calculated from a polynomial 
approximation (MLII) of the form 

Y = 1/A (.5001 + .1649A - .0544A^) (14) 


for 0 < A < .577 


" _ 8.899 + 9.060A + .977A^ (15) 

' ■ A(17.80 + 11.97A + AM 


for .577 < A < 17 


As illustrated in Figure 7 , there is no significant dif- 
ference between the two ML methods for A < 1. The calculated 
values of parameter A listed in Table 2 are all less than 0.4; 
therefore, either ML method would be appropriate for estimation 
of Y from the sample data used in this study. For large sample 
size (m >> 20) and Y < 4, which characterizes this data sample, 
the ML method is favored over the moments method (Pef . 9 ) . 

Comparisons of observed and theoretical distributions of 
gust and gust length are illustrated in Figures 8 through 31. 
The distributions based on ML have smaller values of the shape 
parameter, y; these distributions are more skewed than those 
based on M statistics and deviate more from the observed 
distributions at the extreme percentiles. 

2. Variability of Gynma Distribution Parameters . Vari- 
ability of the estimates of parameters y and B is an indication 
of the variability of the theoretical gust distribution. As 
indicated earlier, y determines the shape of the distribution 
function and 8 is a scaling parameter. Gust percentiles are 
inversely related to $, or directly related to $* , where 
B* - 1/B. 


The variability of y and 6* as a function of filter cut- 
off wavelength, y , and altitude is illustrated in Figures 12 

W 

and 13 for u component gust and in Figures 14 and 15 for v com- 
ponent gust. The parameters were estimated by the moments 
method; the variability is similar for estimates based on the 
method of maximum likelihood. As illustrated in Figures 12 
and 14, the value of y is usually between 2.25 and 3.70 for 
both components; the variability within that range is not 
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Figure 7. Maximum Likelihood Estimates of y Based on Thom (MLI, Ref. 7 ) 
Bury (MLII , Ref. 8). Parameter, A, is Calculated According to 













Table 2. Parameter, A, Calculated from Filtered Jimsphere Data during 
February at Cape Kennedy, Florida 
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♦Insufficient data 
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Figure 10. Observed (Plotted Points) and Theoretical Distributions 
(Curves I, II, and III) of u Component Gust Length, LU, 
at 12 km During February at Cape Kennedy, Florida, for 
= 2470 m 
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Figure 11. Observed (Plotted Points) and Theoretical Distributio 
(Curves I, II, and III) of v Component Gust Length, L 
at 12 km During February at Cape Kennedy, Florida, fo 
X_ = 2470 m 






3* (m/s) 


l.OH 


.8H 


.6H 



100 1,000 10,000 

x^lm) 




Figure 13. Parameter B* of the Gamma Distribution of Zonal Component 
Gust During February at Cape Kennedy as a Function of 
Filter Cut-off Frequency, X , and Altitude 
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Figure 14. Parauneter y of the Ganuna Distribution of Meridional 
Component Gust During February at Cape Kennedy as a 
Function of Filter Cut-off Frequency, X^, and Altitude 





clearly systematic with respect to either altitude or filter 
cut-off frequency. As illustrated in Figures 13 and 15, the 
scaling parameter, 3*, is strongly influenced by filter cut- 
off frequency and, to a somewhat lesser extent, by altitude; 

3* increases as increases; the increase of 3* with altitude 

is most pronounced between 8 and 10 km. 


D. Gamma Distribution of Gust at Reference Altitudes 


1 . Univariate Gamma for Gust Components and Associated 
Gust Lengths . The univariate gamma is calculated from sample 
estimates of the parameters y and 3. In this section, examples 
are presented which illustrate the following; 

• Comparison of theoretical (gamma) and observed 
distributions for various filter cut-off 
wavelengths 

• Seasonal variation of theoretical distributions 

• Altitude variation of theoretical distributions 

• Comparison of zonal and meridional component 
theoretical distributions. 

The parameters used in the theoretical distributions used 
in the various examples were calculated by the moments method 
(described in Section III.C.l). 

Theoretical and observed distributions at 12 km during 
February are illustrated in Figures 16 through 19; all the 
distributions exhibit a similar variation as a function of 
filter cut-off wavelength, X . No large systematic differ- 

ences between the observed and theoretical distributions are 
noted. 

The variation of theoretical distributions with season 
is illustrated in Figures 20 through 23. Three months, 
February, July, and April, were selected to represent the 
winter, summer, and transition seasons, respectively, at Cape 
Kennedy. It is indicated that large seasonal differences for 
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Figure 17. Observed and Theoretical (Ganuna) Distribution of Absolute v Component 
Gust at a Reference Height of 12 km During February at Cape Kennedy 
as a Function of Filter Cut-off Wavelength, X 
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Figure 19. Observed and Theoretical (Gamma) Distribution of v Component Gust 
Length, Lv, at a Reference Height of 12 km During February at Cape 
Kennedy as a Function of Filter Cut-off Wavelength, X 
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Figure 20. Theoretical Distribution of Absolute u Component Gust at 
Reference Height of 12 km During February, April, and 
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Figure 22. Theoretical Distribution of u Component Gust Length 
Lu, at a Reference Height of 12 km During February, 
April, and July 




either u or v component absolute gust are not affected by the 
choice of filter cut-off wavelength. The small seasonal vari- 
ation of gust length illustrated in Figures 22 and 23 is 
probably attributable to the unavoidable inaccuracies in the 
estimation of distribution parameters from sample statistics; 
therefore, there is no significant variation of gust length 
distributions with respect to season. 

The variation of the theoretical distributions of u and v 
component absolute gust with altitude is illustrated in 
Figures 24 and 25. For each filter cut-off frequency, two 
reference altitudes were selected to illustrate the maximum 
variation of the distributions in the 4 to 14 km altitude 
range. The maximum variation does not necessarily occur be- 
tween 4 and ]4 km. It is indicated that there is a significant 
variation of absolute u and v component gust with altitude. 

A comparison of u and v component absolute gust at 8, 12, 
and 14 km during February is illustrated in Figures 26 through 
28. It is indicated that v component gust is consistently 
larger than u component gust for most of the filter cut-off 
wavelengths and altitudes considered. 

2. Gust Component Percentiles . For engineering and design 
applications, it would be desirable to derive an empirical func- 
tion which accurately describes the variation of particular per- 
centiles as a function of filter cut-off wavelength, A^. 

Such a function should be bounded at 0 for A =0 and should 

c 

approach an asymptotic value for large values of A . A func- 
. c 

tion which has this behavior is of the form: 


|X' 


A 

c 

aA + b 
c 


(16) 


where Ix'l^, the gust component percentile, is asymptotic to 

the quantity 1/a and the parameters a and b can be estimated 
by a least-squares technique utilizing a transformation of 
variables in Equation (16) to obtain a linear equation of the 
form 


T = a t bw 


(17) 
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FEB 








Filter Cut-off Wavelengths 



where T and w are the inverse of 


IX'lp and 


respectively. 


A list of parameters a and b, calculated in this manner from 
theoretical percentiles, is given in Table 3 . Curves repre- 
senting Equation 16 for the 50, 90, 95, and 99 percentile gust 
components during February at 12 km over Cape Kennedy are 
illustrated in Figures 29 through 31; for comparison, the 
theoretical percentiles are also shown as plotted symbols. 

It is indicated that Equation 16 provides a very good fit 
to the theoretical percentiles. 


Another function which exhibits the desired asymptotic 
behavior is of the form 


X' 


d e 


-k/Ac 


(18) 


The constants d and k are also estimated by a least- 
squares technique. This function does not fit the observed 
data as well as Equation 16 . 

3. Conditional Gamma Distributions . Conditional pro- 
bability distributions of component gust given gust length for 
u and V component gusts have been calculated by integration of 
the conditional gamma probability density function. Condi- 
tional probabilities of gust, given gust lengths of 0, 100, 

200, 400, and 800 m for A^ = 2470 m, are illustrated in 

Figures 32 and 33 for u and v component gusts, respectively. 

The distribution for zero gust length is a consequence of 
the theoretical model and cannot be supported by aval] able 
data. 

4 . D istribution of Differences Between the Altitudes 
of Zonal and Meridional Gust Components . It follows from the 
definition of gust used in this study that zonal and meridi- 
onal gusts do not necessarily occur at the same altitude. 

The altitude difference, AH, between the gust components is a 
measure of phase difference. As illustrated in Figure 34, 
the distribution of AH can be accurately represented by a 
normal distribution; it is indicated that AH increases as 
filter cut-off wavelength, A , increases. Theoretical dis- 

W 

tributions of AH for February and July at a reference alti- 
tude of 12 km are illustrated in Figure 35; AH is generally 
larger during July especially for larger values of A . The 

variation of AH with altitude is illustrated in Figure 36; 
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Table 


3 


a (s/m) 

1/a (m/s) 
b(s) 


a 

1/a 

b 


a 

1/a 

b 


a 

1/a 

b 


1/a 

b 


a 

1/a 

b 


Parameters a and b of Equation 17 for Calculation 
of Gust Component Percentiles (m/s) at 12 1cm 
During February, April, and July at Cape Kennedy, 
Florida 


50 90 95 99 Percentile 


.0777 

12.87 

909.57 


.1228 

8.14 

592.51 


.0725 

13.80 

1023.6 


.1171 

8.54 

817.06 


.1557 

6.42 

1522.9 


.1358 

7.36 

1510.9 


February 

u Component 

.0461 .0398 

21.68 25.12 

409.57 339.04 

V Component 

.0647 .0552 

15.46 18.12 

276.04 229.65 

April 


u Component 

.0517 .0467 

19.36 21.40 

491.9 411.5 

V Component 

.0681 .0588 

14.69 17.02 

359.43 296.06 

July 

u Component 

.1048 .0927 

9.54 10.79 

705.51 586.51 

V Component 

.0696 .0578 

14.36 17.29 

761.95 648.64 


.0310 

32.20 

244.64 


.0419 

23.87 

167.86 


.0380 

26.35 

304.1 


.0456 

21.93 

212.57 


.0748 

13.38 

427.89 


.0451 

22.16 

484.47 
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Figure 31. Theoretical 95# 95, 90, and 50 Percentile 

Component Gust (Curves Represent Equation 16) 
at 12 km During July at Cape Kennedy, Florida 
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Figure 33. Theoretical Conditional Probability Distribution of v Component 
Gust Given Gust Length, Lv 


Observed Distribution 
Symbol 
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Figure 34. Observed and Theoretical (Normal) Probability Distribution of the Altitude 
Difference, AH, Between u and v Component Gust During February at Cape 
Kennedy (Reference Height = 12 km) 



0009 = Y ^ ^0002 
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Figure 36. Theoretical (Normal) Probability Distribution of the Altitude Difference, AH 
Between u and v Component Gust During February at Cape Kennedy as a Function 
of Reference Altitude 



the theoretical distribution at the various altitudes were 
selected to illustrate the maximum variation of AH distribu- 
tions within the 4 to 14 km altitu'^e range. The variation 
of AH with altitude is not consistent with respect to the 
various filter cut-off wavelengths; i.e., to produce the 
maximum variation with respect to the six reference alti- 
tudes, different altitude combinations are required for each 
type of filtered data. 


E. Distribution of Extreme gust modulus 


Gust, as represented in this analysis, is a two-dimensional 
vector with components, u* and v', and modulus, R, given by 


R = /( u • ) * + (v') ^ 


(19) 


The largest value of R in a gust profile is the extreme gust 
modulus, R__„. A nuir±)er of samples of R„,„, each of size, n, 

rridX 

equal to 150 have been generated from filtered wind profiles 
for the months of February, April, and July. The analysis 
described in this section is concerned with utilizing these 
samples to establish a theoretical probability distribution for 

max. 


It is proposed that the distribution of R is the first 

max 

asymptotic extreme value distribution described by Gumbel (10) ; 
this distribution, which is also known as the Fisher-Tippette 
Type I probability distribution is given by 


•MX) 



( 20 ) 


where y is the reduced variate. 


y = a(X - w) 


( 21 ) 
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The parameters a and p are estimated from the sample mean, 
X, and the sample standard deviation, S . 


a = 




( 22 ) 


where o.. and v„ are population mean and standard deviation which 
N N 

are a function of sample size; for n = 150 Gumbel's tabulation 
(Ref. 10) has = 1.22534 and = .56461. 


The variable X in Equation (20), which herein represents 

R , can be calculated from the reduced variate according to 
max 

the equation 


X 


A- 

- V 


+ 21 
a 


(23) 


Table ( 4 ) contains the values of p and 1/a that have 
been calculated from samples of R during February, April, 

1H3X 

and July as a function of filter wavelength cut-off, 

The function represented by Equation ( 23) appears as a 
straight line when plotted on extreme value probability graph 
paper. The theoretical distributions of during February 

for four filter cut-off wavelengths are illustrated in 
Figure 37; the observed distributions, represented by the 
plotted symbols in Figure 37, do not deviate significantly 
from the theoretical lines. 
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Figure 37. Observed (Plotted Symbols) and Theoretical 
Distribution (Extreme Value) of During 

February at Cape Kennedy for Various Filter 
Cut-off Wavelengths f X 
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SECTION IV. CONCLUSIONS 


] 


The results of this study provide the basis for a vector 
wind model for Cape Kennedy, Florida. A methodology has been 
developed for the derivation and analysis of small-scale per- 
turbations in Jimsphere wind profiles. Gusts in various 
wavelength bands have been derived from these perturbations; 
the probability distribution of gust components and associated 
gust length has been shown to be accurately represented by a 
gamma distribution. Theoretical and observed distributions 
of component gust vary significantly with season, altitude, 
and wavelength range. 
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